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OSCILLATORY GRAVITY WAVES 
IN FLOWING WATER 


Turgut Sarpkaya,! J. M. ASCE 


ABSTRACT 


The main objectives of this study are to determine the conditions of stabil- 
ity of progressive oscillatory gravity waves in water flowing over smooth and 
artifically roughened surfaces and to evaluate the combined effect of various 
interrelated characteristics of flowing water, specifically those of viscous 
shear, turbulent mixing shear, and nonuniform velocity distribution, on their 
propagation. 

At present, there are neither experimental nor conclusive theoretical 
results which constitute a solution to the problem under consideration. It 
thus becomes highly desirable to investigate the conditions of neutral stability 
of progressive oscillatory gravity waves in flowing water and to evaluate the 
effect of various characteristics of flowing water on the velocity of propaga- 
tion of stable waves by a comparison of the experimental results with those 
obtained from the available theories, particularly those of Stokes and Struik. 

The period and various geometrical characteristics of stable ascendant 
oscillatory gravity waves in water flowing in a smooth and an artificially 
roughened rectangular channel were measured for a wide range of depth and 
velocity of flow and suitable settings of wave generator. The following results 
were obtained: 

(1) For a given flow, there is only one neutrally stable wave. All the other 
waves for this flow show either an exponential damping or an amplification 
during their propagation along the channel. It is readily possible to separate 
the two regions by the surface of neutral stability which was determined ex- 
perimentally and presented graphically. 

(2) The celerity of a stable ascendant wave in flowing water is smaller than 
that obtained from the theories of Stokes and Struik. The deviation of the 
measured wave velocity from that obtained from potential theory increases 
with increasing Froude numbers and decreases with increasing Reynolds 
numbers. 

(3) For a channel with an artificially roughened bottom, the deviation is 
greater than the corresponding one in a channel with a smooth bottom. 

(4) The experimentally determined profiles of stable oscillatory gravity 
waves vary about the theoretical profiles as described by Struik. 


INTRODUCTION 
The exact knowledge of the characteristics of waves in flowing water is 
extremely important not only from the scientific point of view but also from 
the practical one. It is particularly important in connection with the 
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penetration of waves into an estuary and probably with the determination of 
scour in the river bed by either ascendant or descendant waves generated by 
wind. The analysis of all these problems requires, first of all, a knowledge 

of the laws governing wave movement in a nonerodible channel. 

Until the last two decades, knowledge of the motion of oscillatory gravity 
waves in water was based largely on theoretical studies. Several rigorous 
solutions, differing chiefly in their results for the motion of the water parti- 
cles involved, were obtained(1, 2, 3, 4). 

Experimental confirmation of these theoretical solutions has been the main 
purpose of recent studies. For the waves in still water, these investiga- 
tions(5, 6, 7) have shown that aside from mass transport the magnitude of the 
difference in value of the wave characteristics defined*by the theories of ro- 
tational and irrotational wave motion is smaller than the experimental error. 
If the waves are generated in a liquid which is already in motion some differ- 
ence would be expected between observations and classical theories resulting 
from the effect of various characteristics of flowing water, specifically those 
of viscous shear, turbulent mixing shear, and nonuniform velocity distribution 
which are definitely interrelated. In the one hand, viscous action in the bound- 
ary layer produces vorticity which diffuses through the fluid and modifies the 
whole velocity distribution. The subsequent wave motion thereby inherits a 
rotational motion, and the usual assumption of a velocity potential becomes 
inadequate. On the other hand, turbulence becomes an important factor and 
it plays a role that may even overshadow that of gravity. Furthermore, both 
the velocity distribution and the turbulence structure of the primary flow are 
completely altered by the wave motion. In addition to these, perhaps one of 
the most complicating factorsis that the waves generated in flowing water are 
not always neutrally stable. Some of them are damped and some of them are 
amplified according to the existing flow conditions. Therefore, the effect of 
the various factors mentioned previously on the propagation of waves in flow- 
ing water can be evaluated, but only after the determination of the conditions 
of neutral stability. 

In this study, an attempt has been made to establish the conditions of stabil- 
ity of ascendant waves in flowing water and to evaluate the effect of the various 
characteristics of flow on their propagation. 


Previous Investigations 


Certain phases of the general problem of the effect of various character- 
istics of the primary flow on the velocity of propagation of disturbances have 
been studied previously. The problem is, however, so complex that none of 
these phases has been completely solved. 

Bazin’s experiments(8) on the propagation of a solitary wave were the first 
to demonstrate the incongruity between observations and the theory despite 
the good agreement in still water. The effect of the nonuniform velocity dis- 
tribution on translation waves has been studied in great detail by Boussines 
Favre(10) has employed the Boussinesq integrals to describe the motion of 
waves occuring during the operation of navigation locks. More than half a cen- 
tury after Boussinesq’s analysis was completed, Keulegan and Patterson pre- 
sented an analysis of the combined effects of turbulent resistance and inertia 
on the motion of flood waves and the other eaves of translation in open chan- 
nels(11). All of these investigations were, it should be remembered, carried 
out for translatory gravity waves in flowing water. 
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Hough(12) presented a solution of certain problems illustrative of the 
effects of viscosity on the motion of the sea, and concluded that the velocity 
of wave propagation is slightly decreased by resistance of the ocean bed, by 
the dissipation of energy through internal viscosity, and by turbulence. 
Abdulla(13) made an attempt towards a solution of wave motion superposed on 
a current which is caused by winds. The current is assumed to follow 
Ekman’s spiral law in speed but to remain constant in direction. From the 
comparison of his results with gravitational waves, Abdullah found that the 
total effect of turbulence and gravity is to give the waves a velocity which is 
smaller than that provided by gravity alone. Biésel(14), has presented a the- 
oretical analysis for the effect of nonuniform velocity distribution on the ve- 
locity of propagation of small amplitude oscillatory gravity waves in open 
channels. This analysis will be discussed later in detail as it partially con- 
firms the experimental results. 

For very long waves, that is for very large values of aL2/d3, in the line- 
arized theory Burns(15) obtained a general equation connecting the wave ve- 
locity with the velocity in the undisturbed stream and solved the equation in 
several special cases. He has shown generally that the velocity of a long 
wave relative to the mean flow is always greater than the value given by the 
classical theory. In the appendix of Burns’ paper, Lighthill examined the 
“critical Froude number” above which no ascendant wave propagation is pos- 
sible. On the basis of some assumptions which will be discussed later, Light- 
hill has found that the nonuniformity in the velocity distribution has opposite 
effect for large and small wave lengths. For Stokian waves where the wave 
motion is predominantly at the surface, the relatively high fluid velocity in 
this region causes the critical Froude number based on the mean velocity to 
be lower. But for really long waves, in which the motion is almost purely 
longitudinal, the decreased velocities near the bottom make a bigger differ- 
ence than do the increased velocities near the surface. 

Considering the second part of the present investigation, one finds that no 
previous theoretical and experimental study has been devoted to the determina- 
tion of the conditions of stability of an oscillatory gravity wave. Because of 
the challenging nature of the stability problem, however, many investigations 
since the time of Boussinesq have been made for the siability of interfacial 
gravity waves, for the stability of parallel flows, and for the stability of dunes 
in the stream beds. Korteweg and de Vries(16) have shown that relatively 
long waves have a tendency to steepen ahead of their peaks, at a rate propor- 
tional toa/d which is resisted by an opposing tendency, independent of ampli- 
tude, and due to the departures from uniformity in the distribution of velocity 
with depth. Korteweg and de Vries have also shown that there is an infinite 
sequence of wave forms, each such that for one particular ratio of a/d to 
d2/L2 the two tendencies to change of form cancel out, leading to a wave of 
permanent type. These waves are cnoidal; more precisely, they are graphs 
of the square of the Jacobian elliptic function cnx for different moduli k1; one 
extreme (kj = 0) is the sinusoidal wave, which occurs when the ratio aL2/d3 
is very small, and at the other extreme (k1 = 1), the wave length becomes 
large and the waves take on more and more the character of a sequence of 
solitary waves. 

For the stability of an elementary wave, starting from a uniform regime, 
the evolution with time of a sinusoidal undulation of very small amplitude was 
studied by H. Jeffreys, and depending upon the case, an exponential damping 
or amplification was found. Recently the stability of two dimensional small 
disturbances in laminar flow on a plane boundary was examined by Yih(17). 
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It was found that for each slope of the channel bottom, a hyperbolic type of 
curve with the arms extending in the directions of wave number 27/L and the 
Reynolds number separates the stable region from the unstable one. 

As has already been mentioned, however, no previous theoretical and 
experimental investigation exists concerning the conditions of stability of an 
oscillatory gravity wave of finite amplitude in flowing water. It thus becomes 
extremely desirable to establish the conditions of stability of waves under 
consideration. 


Theoretical Considerations 


The present status of turbulence theory and the mathematical techniques 
seem to offer little hope that either the effect of turbulence and viscosity or 
the conditions of stability can be directly determined. In general, an oscilla- 
tory gravity wave is unsteady and involves both molar and molecular resist- 
ance. Besides the complexity of wave motion itself, the characteristics of 
turbulence, such as velocity and pressure fluctuations, cannot be asserted a 
priori. However, neglecting the effect of turbulence and considering only the 
effect of viscosity with a nonuniform velocity distribution over the cross sec- 
tion and using the Sommerfeld-Orr equation 18) the problem of stability can 
be formulated conveniently. Despite its simple formulation, the task of evalu- 
ating a rigorous solution of the equations obtained for even a simple velocity 
distribution would be extremely difficult and no attempt has been made at such 
an evaluation. 

For very large Reynolds numbers, the theoretical as well as the practical 
value of this solution seems questionable. As the Reynolds number increases 
beyond the critical range, the eddy viscosity becomes increasingly important. 
In other words, if the Reynolds number is high, the viscous stresses in terms 
of the mean velocity gradients are negligible in comparison with the apparent 
stresses due to momentum transport. 

Despite the difficulties cited above, the effect of nonuniform velocity dis- 
tribution alone on the propagation of stable gravity waves of small amplitude 
can be obtained analytically. In accordance with the above considerations, 
Biesel(14) neglected the viscous terms in Sommerfeld-Orr equation and made 
the wave motion steady by fixing the axes on the wave profile. He further as- 
sumed a linear velocity distribution given by the expression, 


C=G(i+ky) 


and obtained the velocity of an ascendant wave to be 


2 


in which U, and Uj represent the absolute surface and bottom velocities of 
flow, respectively. 

Were this velocity obtained from the superposition of a small-amplitude 
wave velocity on the mean velocity of flow, one would obtain 


Vx 
as m 2 
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If Vt and Vm represent respectively the theoretical and measured wave 
velocities, the deviation ¢ can be defined by the relationship 


4 


m 


If Eq. (2) is used for measured velocities and Eq. (3) for theoretical velocities, 
the deviation ¢ can be obtained as 


Vor 


The results obtained from this equation and their comparison with the experi- 
mental results will be discussed later in detail as it partially confirms the 
experimental results. 

Next we turn our attention to the breaking of waves in flowing water. One 
of the most important characteristic features of oscillatory gravity waves in 
flowing water is the continuous deformation of wave shape and the subsequent 
breaking of waves under the influence of somewhat different factors than those 
encountered in still water. The limiting condition beyond which no wave pro- 
pagation exists is determined primarily from the consideration that the wave 
becomes unstable and breaks when the forward velocity of the water particles 
at the crest exceeds the velocity of the wave. In still water this in general 
may be caused by the shoaling of the depth accompanied by the complicated 
effects of backwash from previously broken waves or by yet little understood 
action of an opposing wind. It has been known theoretically for quite some 
time that the asymptotic limit for waves in deep water is a@/L = 0.14. In fairly 
shallow waters, it has been pointed out already by various authors that a good 
approximation is given by the solitary wave theory despite the fact that it is 
not confirmed by the experiments conducted at M. I. T., in sloping channels. 
For intermediate conditions in medium depth, the limiting wave steepness is 
given by Miche(4) as 


2nd 
0.14 Tanh 6 


Experimental results reported by Danel(19) come very close to this theoreti- 
cal limiting curve. 

In flowing water, however, in addition to the effect of the slope of the chan- 
nel bottom which does not cause any shoaling, the retardation of the advance 
of the forward face of the wave by the opposing current and consequent over- 
taking of the forward face by the remainder of the wave becomes a predomi- 
nant factor. Evidently, the mean velocity of flow as well as the particular 
velocity distribution plays an important role in the determination of limiting 
condition. Therefore, the limiting wave steepness is further related to the 
Froude and Reynolds numbers of the primary flow in addition to the ratio d/L. 

The consideration of all these parameters, however, in the theoretical 
determination of the limiting condition is as difficult as the evaluation of the 
stability conditions. It is essential to make simplifying assumptions in order 
to get some insight into the phenomenon. Lighthill(15) made the assumption 
that a solution exists with V = 0, corresponding to a stationary wave, above 
which no upstream propagation is possible. The consideration of the equation 
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of motion for a stationary wave of sufficiently small amplitude with the par- 
ticular velocity distribution 


and with the necessary boundary conditions leads to the solution for the 


critical Froude number as follows: 
¢ V9 2nd 2 d 
Md (i+é) (4f ) 


To show the effect of nonuniformity in velocity distribution, if the ratio of 
this expression to the classical expression (its limit as€é~0) for uniform flow 
is taken, one finds 


1 1 | 


The fundamental assumption made by Lighthill in the derivation of the fore- 
going equations is not, however, fully acceptable. For the wave propagation 
becomes impossible with the complete breaking of the wave when the velocity 
of flow becomes equal to the celerity of energy transmission, which is in gen- 
eral smaller than the wave velocity except for the case of anomalous disper- 
sion. For the waves of small amplitude, the wave energy is transmitted by 


the group velocity Cg 
4nd/v L 2 ad 


The validity of the foregoing is limited, strictly speaking, to the waves of suf- 
ficiently small amplitude. The assumption is made that the energy transmis- 
. sion for the waves under consideration is of the same nature as that of the 
small amplitude theory. Although, the effect of Froude and Reynolds numbers 
can be included with the help of the deviation 2, for the sake of simplicity, 
this will be omitted. 
If the assumption is made that no wave propagation is possible when the 
celerity of energy transmission is equal to the mean velocity of flow, the 
critical Froude number becomes 


Amd 


It can easily be shown that the condition of the equality of the mean velocity 
of flow to the celerity of energy transmission corresponds precisely to the 
equality of the wave lengths of the two ascendant waves of the same period 
propagating in the same flow. Here it should be noticed that, in general there 
are two different waves of the same period propagating in the same flow. 


Dimensional Analysis 


The basic variables pertinent to this phenomenon are: Amplitude, length 
and velocity of a stable wave; dynamic viscosity and density of the fluid; depth 
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and mean velocity of flow; shape factor of the flume; spacing and the dimen- 
sions of roughness elements; and the acceleration of gravity. These variables 
yield the following functional relationships: 


Od 2,8} 12 


To make the comparison of the wave velocities obtained experimentally with 

those obtained from the theory of irrotational motion, it becomes necessary > 
to introduce the theoretical wave velocity V¢ into the analysis. For an ascend- 

ant wave, Vt is given by the relation 


If both sides of this equation are divided by V,, and the numerator and the 
4 denominator of the right side by gd also, one obtains 
vad 14 
Vin 


Ym 


U vad 


Since C/vgd is a function of only @/d and L/d and in turn a function of only 
the independent dimensionless parameters, V,,/U in Eq. (12) can be replaced 
by Vt/Vm or by the deviation~. If one replaces a/d by a/L and L/d by d/L 
and also d/k by k/d for convenience, the relationships take the form 


a 
L 
da = U Ud 15 


The dependent dimensionless parameters on the left side of the equations 
represent respectively the wave steepness, the ratio d/L, and the deviation. 
The parameters on the right side of the equations represent, in the order of 
their appearance, the Froude number, the Reynolds number, the relative 
roughness, the relative spacing of the roughness elements, and a channel 
shape factor. 

The consideration of all these variables in the determination of the stability 
conditions and in the evaluation of the deviation { is almost impossible. A 
useful result can be attained only in isolated instances and then only if all but 
the most essential parameters are ignored. The existing experimental data 
and the theoretical knowledge are not sufficient to determine the characteris- 
tics of the shape factor 8. Furthermore, if the sides and bottom of the chan- 
nel are smooth, d/k and A/k need not be considered. Thus, for smooth chan- 
nels of rectangular cross section, the relationships take the forms, 
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For the case in which the channel bottom is artificially roughened, 


a 
i =H {FoR Lop} 


In the present investigation, only one type of strip roughness is used; i. e., 
A/k is constant. Therefore: 


a 

=, 18 


The experimental results will be evaluated according to Eqs. (16) and (18). 


Apparatus and Procedure 


Wave Flume 


The glass-walled flume used in these experiments has a rectangular work- 
ing section 85 feet long, 2 feet wide, and 1.5 feet deep. Depth and velocity 
measuring devices and probes for the wave recorder with necessary electri- 
cal circuits were mounted on two aluminum carriages which run on two paral- 
lel steel rails attached to the top of the walls of the flume. 

Artificial roughness consisted of 1/4" x 1/4" continuous brass strips, 
placed at 1" intervals on the steel bottom of the flume perpendicular to the 
longitudinal axis of the flume. The relative spacing )/k of the strips was 4. 
The depth of the flow was measured to the middle of the height of the strips, 

i. e., 0.011' was subtracted from the depth measured to the smooth bottom. 
Actually, the true depth would be somewhere between this depth and the depth 
to the top of the brass strips. 

Water was supplied from the circulation system of the laboratory, and the 
discharge was measured by means of a calibrated venturi meter. in addition, 
a triangular weir, which can be conveniently placed at the end of the flume, 
was used to measure small discharges. 


Wave Generator 


It was necessary, to design a new wave generator which must at once allow 
the passage of flow, produce as nearly as possible the water particle motion 
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in accordance with the theory under consideration, cause a minimum amount 
of additional turbulence, and be simply constructed and easily operated. 

The wave generator developed in the Iowa Institute of Hydraulic Research 
and used throughout the experiments is principally a combination of Rans- 
ford’s classical three-bar linkage or ballistic pendulum with a system similar 
to a venetian blind which is made of metal slats that can be set together at any 
angle to regulate the passage of water only in one direction (Fig. 1). The use 
of metal slats to allow the passage of flow was an idea originated at the 
St. Anthony Falls Hydraulic Laboratory. The venetian-blind type of wave 
generator adequately produces in varying depths of flow the waves under in- 
vestigation throughout the complete range of study and is, moreover, very 
simple to operate. 

In order to eliminate wave reflections from the end of the flume opposite 
the generator, it was necessary to install wave dampers across the section. 

A basket made of No. 32 wire screen was filled with 1/4 and 1/2 inch gravel. 
The basket was hung on two parallel steel rails at a slope of approximately 

30 in the direction of flow. The elevation was so adjusted that the bottom of 
the screen basket was a half wave height below the mean elevation of the water 
surface. This, with one or two wire screens stretched below the section, 
provided excellent wave absorption. 


Experimental Procedure 


The amplitude, length, and shape of waves were determined from photo- 
graphs. Readings were made to thousandths of a foot under the microscope 


with no difficulty. 
Wave periods were obtained from recordings made with a Sanborn Twin- 


Viso recorder. Wave velocities were also obtained by placing two parallel 
wire gages 10 feet apart in the wave tank and recording the wave profile from 
these gages on the Sanborn Twin-Viso recorder. A timing device made four 
marks per second on the recorded permapaper. Thus, knowing the distance 
between the gages and the time required for the wave to travel from one gage 
to the next, the wave velocity between the gages was determined. 

Velocity profiles on the axis of the flume were obtained with a 1/8" diam- 
eter pitot tube connected to two manometers which were set at a slope of 45° 
from the vertical, in order to magnify the differences V2 times. 

The general experimental procedure in the tests described herein was as 
follows: The discharge was held constant for varying depths of flow. For 
each depth the setting of the wave generator was varied several times, and 
wave heights were recorded at different sections of the flume with the same 
wave recorder and compared. Necessary adjustments of the generator, such 
as period and the positions of the connecting rods, were made until a stable 
wave was obtained. After the stability of wave had been established the char- 
acteristics of the resulting waves were photographed and recorded. 


Discussion of Results 


Stability of Waves 


An oscillatory gravity wave is said to be stable when the matter, energy 
and the phase transported remain constant throughout the propagation of wave. 
Therefore, instability corresponds to the continuous variation of the foregoing 


parameters. 
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It has already been pointed out by various authors that the assumption of 
constant power transmission is not quite admissible particularly for shallow 
areas of flat slope and, reductions of as much as 30% in wave heights may 
occur due to bottom friction in still water. When the media of wave propaga- 
tion is already in motion, however, additional energy dissipating factors such 
as internal viscous shear and turbulent mixing shear come into consideration. 
The waves can remain stable then if only the energy withdrawn from the waves 
by means of these energy dissipating agencies is exactly balanced with the 
energy added by the flow. Evidently the propagation of waves alters com- 
pletely both the velocity distribution and the turbulence structure of the pri- 
mary flow. The supplied and the dissipated energies which are equal origi- 
nally, become unbalanced because of the complete alteration of the character - 
istics of the primary flow and a certain amount of energy, the difference be- 
tween the supplied and the dissipated energies, is left over. This unbalanced 
energy is used, in general, either to amplify or to damp the waves. There- 
fore, from the energy transfer point of view, instability corresponds to an ex- 
cess or deficiency of the supplied energy with respect to the dissipated energy 
by. various resistances. 

The parameters which change the characteristics of the primary flow are 
essentially the oscillations of the laminar boundary layer due to periodic vari- 
ations of the pressure, the orbital motion of the water particles due to wave 
movement and the mass transport and, as a consequence of the latter, a mo- 
mentum transport in the direction of wave motion. The momentum trans- 
ported per wave length is simply equal to wave velocity times the mass trans- 
ported during one wave period. Recent studies have shown that the disturb- 
ances on the water surface even of small amplitude, change the mixing length 
of turbulence considerably. It is an empirically known fact that the turbulence 
creates a resistive force roughly dependent upon some power of the velocity 
and indirectly upon some power of the hydraulic radius with a proportionality 
constant depending upon the stream bed. This empirical formula could per- 
haps be more rationalized after the determination of the variation of the inten- 
sity and scale of turbulence as a function of Froude and Reynolds numbers, so 
that the wave energy scattered by turbulence can be formulated. Subsequently 
the attenuation of wave which is at least of the same order of magnitude as 
the energy scattered can be evaluated. 

Because the determination of the stability conditions lacks a sound theoreti- 
cal solution, the direct approach by way of experiment appears to be the unique 
method, and may be expected to yield satisfactory results. In order to obtain 
the conditions of stability of oscillatory gravity waves in water flowing in a 
smooth channel, the steepness and the ratio d/L of stable waves were deter- 
mined experimentally and plotted versus. the Reynolds number of the primary 
flow. A series of curves was obtained for constant values of Froude number 
as is shown in Figs. 2 and 3. Each curve is similar in shape to the one ob- 
tained by Yih for the small-amplitude waves in laminar flow. In Yih’s solu- 
tion for each slope of the channel bottom, a hyperbolic type of curve with the 
arms extending in the directions of the wave number and the Reynolds number 
separates the stable region trom the unstable one. From Figs. 2 and 3, two 
new graphs were obtained one showing the curves of constant Reynolds num- 
bers for L/d versus Froude number (Fig. 4) and the other showing the sur- 
face of neutral stability as a function of L/d, Froude number, and the Reynolds 
number (Fig. 5). Figures 2, 3, 4, and 5 show clearly that for a pair of given 
values of the Froude and Reynolds numbers; i. e., for a given flow, there is 
only one wave which is neutrally stable. For this flow, the other waves are 
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either damped (below the curves of constant Froude numbers) or amplified 
(above the curves of constant Froude numbers) as they propagate along the 
channel. A wave which is amplified along its course cannot, obviously grow 
indefinitely. After a certain distance throughout which the wave grows con- 
tinuously in amplitude, breaking occurs near the crest, accompanied by the 
appearance of foam extending down the front slope of the wave. The foam 
appears to be detached from the mass of the water and rolls over the advanc- 
ing face. The subsequent development of the wave was not examined. Thus, 
the possibility of stabilization of a wave by breaking is excluded from the 
consideration herein. 

For a given flow, the ratio d/L and the steepness of the stable wave which 
would propagate on this current can easily be approximated from Figs. 2 
and 3. Figure 2 also reveals that for a given Froude number as the amplitude 
of the wave increases, the critical Reynolds number for the stable wave de- 
creases. But for large values of the Reynolds number, the steepness of the 
stable wave depends only on the Froude number. As can be seen clearly from 
Fig 3, the ratio d/L depends on both the Reynolds and Froude numbers. How- 
ever, the dependency of the ratio d/L on the Reynolds number diminishes 
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with increasing Reynolds numbers if the Froude number is kept constant. 


Evaluation of the Deviation in Wave Velocity 


For a preliminary comparison of the results, Laplace-Airy equation can 
conveniently be taken into consideration: 


ge 2nd 
C= 2n L 


If the mean velocity of flow is added to the wave velocity, the celerity of an 
ascendant wave is obtained: 


If Eqs. (19) and (20) are divided by vgd, one obtains 


22 


and 


The computed values of wave celerity from measured values of wave length, 
wave period, and mean velocity of flow, for each wave studied, are compared 
with the corresponding values computed from Eq. (21). The data, presented 
graphically in Fig. 6, show that ascendant waves propagate consistently more 
slowly in flowing water than the Laplace-Airy theory would indicate. The sig- 
nificance of the comparison is, however, limited by the fact that in the 
Laplace-Airy theory the effect of wave steepness on the celerity of wave is 
not taken into account. 

In the third approximation, the celerity of a wave Ct, according to Struik, 
is given by the following equation, 


ct= {pb tanh 249} 23 


in which 


Thus the velocity is a function of the wave steepness as well as of the ratio 
d/L. The factor A becomes increasingly important as the ratio d/L becomes 
smaller. In order to show this importance and to facilitate the computation 
of Ct, a parametera, which is defined by the expression 


25 


Cano 
is introduced into the analysis: The values of a are taken from Fig. 7 for 
corresponding values of wave steepness and the ratio d/L , and the values of 
Cq = 0 are computed from Eq. (19). For the same parameters, Fig. 8 gives 
the values of Ct directly, although less reliably because of the large variation 
of the celerity with wave steepness. 


815-16 


+ 


seo 
oso 


| | | | — 
| 
\ 
815-17 


a 
0900 | = 
\ 
—— 
0v00 
| | Meg 
| | | 
| | | \ | 
| | \ 
020°0 | | | | 
| 4 
| | | 
| | \ 
| 
| | 
| | | | \ 
| | 
96000 | 
| | | E\ 
92600. | 
| 
| | \ 
| | | | 
| 
815-18 


Brecking Limit 
(Miche) 


|_| 


+> 


| 


| 
0.15 | 
| | | 
dit A 
| j 
| | | 
| | | | ° | 
] | 
T 
| 
| Ea 0.075 
—0.045 
0.035 
| 
| 
0.03 | | | | 
| | | | — 
©0035 | 
t | 20.04 
@ 0.088 
©0065 
| | ® 0.075 | 
AS A PUNCTION OF a/L and d 
815-19 


YU NOILONA V SVQ 6 


@ 
4 


| 


| 
| | | | 
| 660000 
3 
| | | | 
| 
| | | 3 | | | 
| | 
| | | | 
| lo 
| ~ 
| | } 
| | | 
| , le | | 
| | | 
| 
| 
| | | } 
| iy | | 
| | | 
| @ | 
| | 
1 o 
| 
| j | 
| | | | | | 
0/0 o | | | | | 
815-20 


After the fundamental theoretical solution had been properly prepared 
for comparison with the experimental results, the aforementioned { which 
was defined by the equation 


was computed and plotted versus Reynolds number and a series of curves 
were obtained for constant Froude numbers (Fig. 9). As found from dimen- 
sional analysis, the ¢ is a function of the Froude and Reynolds numbers. 
Figure 9 shows that for a fixed Froude number, the dependency of the devia- 
tion on the Reynolds number diminishes as the Reynolds number increases. 

In order to compare the values of { determined from the experiments with 
those obtained from Biésel’s analysis for sinusoidal waves, the ratio d/L was 
plotted versus the Froude number for neutrally stable waves, and a series of 
curves of constant £ was obtained in Fig. 10. Equation (5) is presented in the 
same manner in Fig. 11, for K = 1.34, which was determined from the velo- 
city profiles obtained experimentally. The similarity in the trends of the ex- 
perimental and the theoretical results is obvious. The agreement is not re- 
markable. This, however, is not unexpected, since in the analysis turbulence 
was not taken into consideration, and the waves were of infinitesimal ampli- 
tude, and as a consequence of the latter limitation, the mass transport in the 
direction of wave motion was out of question. 

In Fig. 8 the deviation ¢ is presented as a function of the ratio d/L and 
the wave steepness. For the values of d/L less than approximately 0.40, ¢ 
increases rapidly with increasing wave steepness. As the ratio d/L becomes 
larger, the wave steepness becomes the dominant factor. Here it should be 
remembered that the foregoing comparison and conclusions have so far been 
made only for experiments in which the channel walls and bottom were smooth. 
For a channel with an artificially roughened bottom and smooth sides, similar 
plots are presented in Figs. 12, 13, and 14. These figures clearly show that 
the deviation for the roughened channel is greater than the corresponding one 
for the smooth channel. Unfortunately, it is not possible to draw a definite 
conclusion about the dependency of the deviation on the Reynolds number. 
Moreover, a careful examination of Figs. 14 and 15 reveals that for a given 
Froude number the deviation does not vary greatly with the degree of rough- 
ness. Therefore, the dependency of the deviation on the Reynolds number is 
slight. However, the question remains open as to whether, in the case of an 
artificially roughened channel, the deviation is a function of the Froude num- 
ber alone or depends also on the Reynolds number. This point could be clari- 
fied through experiments with new roughness elements geometrically similar 
to the one previously used but of different size. Such experiments would en- 
able one to prepare a family of plots with a constant Froude number for each 
plot and each one containing a series of curves for constant relative rough- 
nesses as a function of the deviation and the Reynolds number. Figure 15 
presents such a family of curves, although incomplete, for a Froude number 
of 0.30. 


Relationship between Ascendant and Descendant Waves 


If waves are generated in a fluid which is already in motion by placing the 
wave generator somewhere near the middle of the flume, one observes that 
descendant waves are much longer than ascendant ones whereas their period 


is equal. 
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According to the Laplace-Airy solution, the ratio of the lengths of an 
ascendant and a descendant wave of the same period can be written as 


follows: 


Constant Froude numbers are plotted versus Lg/d and L,/d in Fig. 16. This 
figure shows that a descendant wave propagates faster than an ascendant wave 
provided that they both have the same period and travel in the same flow. 
This fact becomes more apparent as the Froude number is increased and 
shows the difficulty of testing descendant waves in a short test channel. The 
limiting barrier beyond which no ascendant wave propagation is possible 

(Eq. 11) is also shown in the same figure. 


Comparison of Theoretical and Experimental Wave Profiles 


Typical experimental profiles selected from the numerous runs are com- 
pared with the theoretical profiles given by the equation (3): 


+m (1+ 2 cothind)}8 cos mx 


coth md )8 cos 2mx 


+a “(3 +1) 8 cos3mx 26 


in which § is determined from the equation: 
a=28{i+ cothimd+i3 coth* cot#tmd—)} 27 


Although no theoretical function represents the wave profile exactly for all 
possible combinations of ratio d/L and wave steepness, it was found that the 
above equation gives the most consistent agreement between observations and 
theory within the range tested. The deviation of the experimental profiles 
from the theory was as follows: For large amplitudes, the experimental pro- 
files have a sharper crest with more gentle side slopes and a flatter trough 
than given by the theory, while the lower amplitudes exhibit characteristics 
of an opposite nature. 

Thus, for the amplitudes and depths examined in this investigation, the 
experimental profiles vary about the theoretical profile as described by 
Struik. The wave profile in Fig. 17 is representative of the numerous profiles 


tested. 


CONCLUSIONS 


The detailed laboratory and theoretical investigation of oscillatory gravity 
waves in flowing water described herein produced results from which the 
following conclusions can be drawn: 
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(1) For a given flow, there is only one neutrally stable wave. All the other 
waves for this flow show either an exponential damping or an amplification 
during their propagation along the channel. It is readily possible to separate 
the two regions by the surface of neutral stability which was determined ex- 
perimentally and presented graphically. In view of the complexity of stability 
phenomena, it is no wonder that there is considerable difficulty in finding a 
reasonable solution for the stability of oscillatory gravity waves in flowing 
water. This author has no basic solution to present and leaves the theoretical 
determination of this elusive problem to future investigators. 

(2) The celerity of a stable ascendant wave in flowing water is up to 10% 
smaller than that obtained from the theories of Stokes and Struik. This find- 
ing is contrary to the commonly held assumption that the characteristics of 
waves propagated in flowing water are identical with those propagated in still 
water. Thus, the rule of superimposition is not permissible. Two present 
theories are found to be inadequate for the waves in flowing water. 

The deviation of the measured wave velocity from that obtained from po- 
tential theory increases with increasing Froude numbers and decreases with 
increasing Reynolds numbers. The larger the Reynolds number is, the small- 
er is the dependency of the deviation on the Reynolds number. 

(3) For a channel with an artificially roughened bottom, the deviation is 
greater than the corresponding one in a channel with a smooth bottom. For 
example, for a Froude number 0.30 and Reynolds number 105, deviation is 
about 6% for a smooth channel and about 9% for a rough channel. For rough 
channels, however, no conclusion can be drawn as to the dependency of the 
deviation on the Reynolds number because of experimental limitations. Indi- 
cations are such that the Reynolds number has either little or no effect on the 
deviation. 

(4) The experimentally determined profiles of stable ascendant oscillatory 
gravity waves vary about the theoretical profiles as described by Struik. It 
was found that Struik’s profile equation gives the most consistent agreement 
between observations and theory within the range tested. 

(5) The hydraulic engineer who finds himself concerned with the problems 
of oscillatory gravity waves in flowing water should consider this paper as 
but a basic introduction to a body of knowledge and scant available literature 
that exists on wave motion in flowing water. There remain multitudes of un- 
solved problems of both a theoretical and a practical nature, but a high degree 
of understanding of the mechanics involved and an appreciation of the impor- 
tance of sound theoretical concepts must be brought to their solution. 
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